
fvVectorMatrix UEqn
(
    fvm::ddt( U )
    + fvm::div( phi, U )
    - fvm::laplacian( nu, U )
);

solve( UEqn == -fvc::grad( p ) );
